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1.  INTRODUCTION 


A.  PURPOSE  AND  MOTIVATION 

The  purpose  ot  this  study  is  to  examine  the  effects  of  solar  radiation  and  its 
attenuation  on  the  dynamics  of  the  ocean  mixed  layer  and  to  improve  the  accuracy  of 
short-wave  radiation  parameterization  schemes  to  be  used  in  mixed  layer  predictions.  In 
this  context,  "attenuation"  is  synonymous  with  the  term  "extinction,”  referring  to  the 
combined  effects  of  absorption  and  scattering  of  electromagnetic  energy. 

The  importance  of  accurate  analyses  and  predictions  of  upper  ocean  structure  in  the 
fields  of  ocean  acoustics  and  global  heat  budget  analysis  is  well  documented.  The  rapid 
evolution  of  micro-  and  mini-computers  ha«  recently  provided  the  capability  to  evaluate 
and  predict  the  upper  (xean  thermal  structure  using  models  and  databases  locally  available 
to  the  operator.  With  direct  access  to  satellite  observations  providing  an  additional  real¬ 
time  data  source  for  these  prediction  systems,  ocean  optical  characteristics  are  better 
known.  Mixed  layer  mtxlels  need,  therefore,  to  be  able  to  parameterize  more  accurately 
the  effect  of  (Kean  optical  characteristics  in  the  dynamic  mixing  processes  of  the  upper 
(Kean. 

One-dimensional  mixed  layer  models,  such  as  those  described  by  Garwood  (1977) 
and  Mellor  and  Yamada  (1977),  are  typical  of  the  models  available  for  on-scene 
prediction  systems.  The  numerical  complexity  of  three-dimensional  models  limits  their 
use  in  the  smaller  computers  available  for  tactical  acoustic  forecasts  performed  kKally. 


Current  one-dimensional  mixed  layer  models  use  simple  parameterization  schemes  to 
prescrilie  the  attenuation  of  solar  radiation  in  the  ocean.  These  parameterizations 
generally  include  simplifying  assumptions  that  may  ignore  some  physical  processes 
important  to  the  evolution  of  the  mixed  layer.  To  take  advantage  of  improved  optical 
data  sources  such  as  colorimetric  satellites,  improvement  of  these  parameterizations  is 
required. 

This  study  examines  various  improvements  to  these  parameterization  schemes  and 
demonstrates  their  effect  on  analytical  and  numerical  representations  of  mixed  layer 
dynamics. 

B.  BACKGROUND 

The  field  of  ocean  optics  has  developed  as  oceanographers  and  mariners  have 
attempted  to  explain  the  wide  variability  of  ocean  color  and  visibility.  Thus,  an  accurate 
and  complex  description  of  the  various  effects  of  the  ocean  on  solar  radiation  exists. 
Jerlov  (1976)  provides  a  good  overview  of  ocean  optics  and  describes  reflection, 
refraction  and  attenuation,  the  primary  interactions  of  sunlight  with  the  ocean  and 
atmosphere. 

Ocean  mixed  layer  modeling  is  a  relatively  new  field  in  oceanography.  The 
majority  of  the  work  in  the  field  has  emphasized  the  development  of  accu’^ate  schemes 
to  describe  the  evolution  of  temperature  and  mixed  layer  depth  as  a  result  of  energy 
exchanges  between  the  atmosphere  and  ocean.  The  absorption  of  solar  radiation  is  known 
to  be  a  major  determinant  of  the  mixed  layer  thermal  structure  and  is  described  by  the 
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vertical  distribution  of  the  downward  solar  irradiance  in  the  upper  ocean. 

The  magnitude  of  upward  irradiance  in  the  ocean,  a  product  of  back-scattering  by 
water  molecules  and  suspended  particulates,  is  known  to  be  only  a  small  fraction  of  the 
downward  irradiance  and  is  commonly  neglected.  The  amount  of  solar  radiation 
scattered  out  of  a  vertical  column  of  water  is  approximately  equal  to  the  amount  scattered 
into  the  column,  thus  the  vertical  attenuation  of  downward  irradiance  is  used  to  represent 
the  absorption  of  solar  radiation  as  heat  energy. 

The  functional  form  most  commonly  used  for  the  attenuation  of  sunlight  in  the 
ocean  is  the  exponential  relationship 

/(z)  =  (II) 

where  z  is  the  vertical  axis  defined  positive  upward  from  the  ocean  surface,  k  is  the 
attenuation  coefficient,  /  is  the  solar  irradiance  at  depth  z.  and  k  is  the  net  solar 
irradiance  just  below  the  surface  (insolation  at  the  surface  minus  the  fraction  lost  due  to 
the  albedo  of  the  (Kean  surface).  Albedo  includes  the  combined  effects  of  surface 
reflectance  and  subsurface  back-scattering  out  of  the  (Kean. 

The  divergence  of  downward  irradiance  (Q)  is  the  amount  of  solar  radiation  that  is 
absorbed  over  depth  and  can  be  written  as 

CKz)  -  —  -  =  kl  (1-2) 

dz 

Different  methods  of  determining  both  /  and  /o  are  analyzed  in  this  study.  The 
optical  characteristics  of  the  atmosphere  are  the  primary  determinant  of  /o,  but  only  the 
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role  of  the  ocean’s  optical  characteristics  are  examined  in  this  study.  It  is  impossible, 
however,  to  totally  separate  the  oceanic  and  atmospheric  optical  properties  in  transiting 
the  air-sea  interface.  Hence,  certain  assumptions  and  generalizations  are  made  regarding 
the  atmospheric  effects  so  that  the  oceanic  effects  can  be  isolated  and  examined  in  detail. 
These  assumptions  will  be  described  as  needed  in  the  following  chapters. 
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II.  OPTICAL  ANALYSIS 


In  this  analysis,  the  penetration  of  solar  radiation  through  the  ocean  surface  is 
examined  in  detail.  Three  distinct  optical  effects,  albedo,  refraction,  and  attenuation 
influence  the  solar  radiation  incident  upon  the  ocean  surface.  Combining  concepts  from 
the  field  of  actinometry  with  ocean  optics,  an  accurate  account  is  made  of  the  oceanic 
processes  that  influence  solar  radiation. 

A.  INSOLATION 

Insolation  is  defined  as  the  amount  of  solar  radiation  per  unit  time  incident  on  a 
horizontal  plane  at  a  given  location  on  the  earth’s  surface.  Many  procedures  have  been 
developed  to  calculate  insolation  at  a  specific  point  on  the  earth’s  surface  based  on  the 
relative  position  of  the  sun  and  on  the  atmospheric  conditions.  Methods  currently  used 
to  incorporate  atmospheric  effects  into  ocean  mixed  layer  studies  are  approximate,  as 
cloud  cover  information  is  frequently  the  only  meteorological  data  recorded. 
Atmospheric  transmittance  (f)  must  be  known  to  accurately  examine  the  atmospheric 
influence.  Atmospheric  transmittance  is  defined 


where  /t  is  the  irradiance  without  an  atmosphere.  It  is  a  function  of  solar  altitude  {<i>)  and 
earth-sun  separation.  The  term  /o,  is  the  irradiance  just  above  the  ocean  surface.  The 
value  of  r  can  be  obtained  for  a  specific  time  and  location  through  measurements  of  /o. 
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with  an  instrument  such  as  a  pyranometer  and  analytical  computation  of  A.  Page  (1986) 
provides  algorithms  used  in  this  study  to  calculate  <i>  based  on  local  time,  date  and 
location.  Atmospheric  soundings  that  provide  measurements  of  water  vapor  and 
temperature  can  be  used  to  estimate  F,  and  models  such  as  LOWTRAN’  incorporate  the 
attenuation  effect  of  aerosols  on  Iq.  Remote  sensing  methods  to  estimate  insolation  from 
satellite  measurements  of  atmospheric  constituents  have  been  developed.  Clifford  and 
Hay  (1984)  describe  models  that  use  satellite  data  to  predict  atmospheric  attenuation. 
Although  improved  methods  to  describe  atmospheric  effects  on  insolation  predictions  need 
to  be  incorporated  in  mixed  layer  studies,  this  study  focuses  only  on  optical  effects  of 
the  ocean  using  previously  developed  methods  to  estimate  insolation. 

The  two  methods  most  commonly  used  to  estimate  insolation  values  in  mixed  layer 
studies  are  the  empirically  derived  formulae  described  by  Reed  (1977)  and  Lumb  (1964). 
Reed  uses  a  formula  developed  by  Seckel  and  Beaudry  (1973),  also  known  as  the 
Smithsonian  formula,  to  calculate  clear  sky  mean  daily  insolation  (/oo).  The  formula, 
with  coefficient  values  for  two  latitude  belts  to  calculate  Ad  in  W/m*  is 


^  Aq cose  +  sine  +  A2CX)s2€  +  B2sin2e 


Latitude  20°$  to  40°N 
Ao  =  -15.82  +  326.87cosL 

A,  -  9.63+192.44cos(L+9(y’) 

B,  =  -3.27+ 108. 70sinL 

A,  =  -0.64  +  7.80sin(2L-90*’) 

B:  =  -0.50+14.42cos(2A-10") 


Latitude  40°N  to  60°N 
Ao  =  342.61-1. 97L-0.018L’ 
A,  =  52.08-5. 86L+0.043L^ 
5,  =  -4.80+2.46L-0.017L’ 

A,  =  1.08-0.47L+0.011L^ 

B,  =  -38.79 +  2.43L-0.034L^ 


'Mtxlel  developed  by  the  U  S.  Air  Force  Geophysics  Laboratory  to  calculate 
absorption  along  atmospheric  paths  as  a  function  of  wavelength. 
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where  e  =  (Julian  dale  -  21  >(360/365)  and  L  is  the  latitude.  Hourly  clear  sky  irradiance 

Uoc)  is 

Joe  =  4i)«n4>  (2-3) 

where  4>  is  the  mean  value  of  the  solar  altitude  for  the  particular  hour.  Negative  <#»  values 
are  converted  to  zero  to  prohibit  negative  insrrlation  values  when  the  sun  is  below  the 
horizon.  Atmospheric  adjustment  based  on  total  sky  cloud  cover  is  provided  by 
Laevastu’s  (1960)  formula: 

/o,  =  1^0  -0.60C^)  (2.4) 

where  C  is  the  fractional  amount  of  cloud  cover  resolved  to  the  nearest  tenth. 

Lumb  (1964)  incorporates  cloud  type  and  vertical  cloud  distribution  (low,  medium 
and  high)  as  well  as  total  sky  cloud  cover: 

fp.  =  S{a  +  5sin4»)sin4>  (2.5) 

where  5  is  the  solar  constant,  representing  the  total  irradiance  at  the  outer  limit  of  the 
earth’s  atmosphere  at  the  mean  earth  -  sun  separation.  Measurements  by  Wallace  and 
Hobbs  (1977)  indicate  S  is  1380  W/m^  as  opposed  to  the  value  of  1350  W/m^  originally 
used  by  Lumb.  The  empirical  coefficients  a  and  b  vary  according  to  the  cloud 
conditions.  For  relatively  clear  skies,  a  =  0.61  and  b=0.20. 

Lumb’s  method  is  generally  credited  with  greater  accuracy  (Simpson  and  Paulson 
1979;  Lind  and  Katsaros,  1986)  if  data  on  cloud  type  as  well  as  cloud  cover  are 
available.  Figure  1  provides  a  comparison  of  the  two  methods  for  different  fractional 
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cloud  cover  values.  Though  the  difference  in  /q  between  the  two  methods  is  considerable, 
and  may  be  the  largest  source  of  error  in  any  calculation  of  ocean  surface  heat  flux,  the 
verification  of  these  methods  is  beyond  the  scope  of  this  study.  The  method  described 
by  Reed  is  used  in  this  study  and  assumed  adequate  for  the  purposes  of  this  study. 

B.  ALBEDO 

Albedo  {R)  is  defined  as  the  ratio  of  upward  irradiance  over  downward  irradiance 
just  above  a  surface  (Jerlov,  1976).  At  the  ocean  surface,  the  upward  irradiance  includes 
both  surface  reflectance  and  emergent  irradiance  from  beneath  the  suriace.  Extensive 
measurements  by  Payne  (1972)  indicated  that  the  emergent  irradiance  ranges  from  0.5  to 
2  percent  of  the  total  albedo.  In  many  atmospheric  and  oceanic  models  that  include 
surface  radiation  flux  considerations,  the  emergent  irradiance  is  ignored,  and  albedo 
consists  only  of  surface  reflectance.  Empirical  methods  for  calculating  albedo  (Payne, 
1972  and  Coakley,  1979)  include  the  emergent  irradiance  component,  because  the 
emergent  component  cannot  be  distinguished  from  the  reflected  component. 

Albedo  is  a  function  of  f,  <t>  and  wind  speed  (V).  The  value  of  T  determines  the 
relative  amount  of  irradiance  from  direct  sunlight  (with  the  sun  considered  as  a  point 
source)  and  from  diffuse  sunlight  or  skylight.  If  the  diffuse  sunlight  is  considered 
isotropic,  then  the  diffuse  albedo  is  not  a  function  of  <f>  and  is  approximately  0.06. 
Albedo  of  the  direct  sunlight  is  a  function  of  0  and  can  be  approximated  by  the  Fresnel 
equation  for  reflectance  of  a  flat  water  surface  for  an  unpolarized  source  (Jerlov,  1976). 
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Figure  1 .  Hourly  surface  irradiance  values  for  different  cloud  cover  conditions  using 
methtxls  described  by  Reed  (1977)  -  solid,  and  Lumb  (1964)  -  dashed.  Curves 
computed  for  Julian  date  171  at  50°n,  145°W. 


^  ^  2  sin^(9  -4>)  tan^(B  -<t>)|  (2.6) 

2  sin^(0  +  4>)  tan^(e  +  <j>)] 

where  r  is  the  reflection  coefficient,  i.e.  the  fractional  amount  of  an  incoming  solar  ray 
that  is  reflected  from  the  surface,  and  6  is  the  refracted  angle  of  a  ray  passing  through 
the  ocean  surface  measured  from  the  horizontal.  The  magnitude  of  6  can  easily  be 
calculated  from  Snell's  law; 
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where  n  is  the  index  of  refraction  of  water  relative  to  the  air  above,  approximately  equal 
to  1 .34  for  seawater. 


For  a  perfectly  transparent  atmosphere  (F  =  1.0),  all  irradiance  is  direct,  and  the 
albedo  varies  strongly  as  a  function  of  4>.  As  F  approaches  0,  the  direct  component 
disappears  and  the  albedo  approaches  0.06.  Thus  the  value  of  F  prescribes  the  sensitivity 
of  the  total  albedo  to  </>.  Some  mixed  layer  experiments,  such  as  the  application  of  the 
Garwcxxl  (1977)  model  to  ocean  station  ”P,"  have  used  a  constant  albedo  value 
throughout  the  day,  ignoring  any  variation  with  4).  Payne  developed  an  empirical  table 
that  provides  albedo  values  for  various  values  of  <l)  and  F.  If  F  can  be  measured  or 
accurately  calculated  for  a  location,  Payne’s  tables  have  proven  to  be  extremely  accurate 
(Simpson  and  Paulson,  1979;  KaLsaros  et  al.,  1985). 

Coakley  (1979)  developed  an  empirical  equation  for  calculating  clear  sky  albedo  as 
a  function  of  </)  as  follows: 


R 


0.05 


(2.8) 


sin4>  +  0.15 

Figure  2  illustrates  the  effect  of  different  albedo  schemes  on  surface  irradiance  for 
a  typical  clear  sky  day  at  ocean  station  "P."  The  scheme  presented  by  Reed  (1977)  was 
used  to  estimate  surface  irradiance  values.  For  this  case,  where  most  of  the  irradiance 
IS  direct,  the  Fresnel  effect  can  be  seen  in  the  Payne  and  Coakley  schemes  as  surface 
reflectance  increases  significantly  for  low  <j)'s,  and  is  much  less  than  .06  at  high 


Since  even  with  a  "clear"  sky,  F  is  less  than  1 .0,  Payne’s  R  values,  which  represent  F 
=  1.0,  may  depend  too  strongly  on  <j>,  and  Coakley’s  formula  should  be  more  accurate 
since  it  is  based  only  on  clear  sky  measurements.  This  hypothesis  is  supported  in  the 
clear  sky  radiation  analyses  of  Behr  (1990).  Coakley’s  formula  is  used  for  the  analysis 


Figure  2.  Clear  sky  surface  irradiance  at  50^,  145®W  on  Julian  date  171  after 
application  of  different  albedo  schemes.  Solid  -  constant  albedo  (.61);  dashed  -  albedo 
by  Payne  (1972);  dotted  -  albedo  by  Coakley  (1979). 


The  effect  of  surface  waves  on  albedo  is  complex,  as  wave  size  and  shape  as  well 
as  white-capping  and  wave  direction  all  influence  the  surface  reflectance.  Studies  by  Cox 
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and  Munk  (1954)  and  Katsaros  et  al.  (1985)  indicate  that  the  effect  of  surface  waves  for 
values  of  0  greater  than  30“  is  negligible  as  the  surface  area  of  wave  faces  tilted  away 
from  the  sun  is  approximately  equal  to  the  surface  area  of  faces  tilted  toward  the  sun. 
At  low  sun  angles,  however,  the  wave  surfaces  tilted  away  from  the  sun  are  shadowed 
by  the  wave  crests,  and  thus  a  reduction  in  effective  angle  between  the  sun  and  the  water 
surface  prevails,  reducing  the  reflectance.  Quantifying  this  reduction  is  difficult,  as  the 
angle  between  wind  direction  and  solar  azimuth,  white-capping  and  wave  shape  are 
important  factors.  Katsaros  et  al.  (1985)  address  this  difficulty  and  offer 
recommendations  for  additional  work  in  the  field.  An  empirically  based  parameterization 
of  the  surface  wave  effect  is  offered  by  Katsaros  et  al.  for  wind  speeds  less  than  12  m/s; 

R  =  (-0.036)  K  (2.9) 

where  /?o  is  the  albedo  without  consideration  of  wind  speed,  and  V  is  the  wind  .speed  in 
m/s. 

Figure  3  illu.strates  the  effect  of  this  parameterization  for  a  typical  winter  day  at 
(Kean  station  "P"  where  low  4>  prevails  throughout  the  day.  Wind  speed  data  were 
extracted  from  10  days  of  meteorological  observations  at  iKean  station  "P"  during  1961 . 
Values  of  V  ranged  from  4  to  10  m/s.  This  wind  speed  parameterization  is  included  in 
the  analyses  to  follow.  Wind  speeds  in  excess  of  12  m/s  are  treated  as  12  m/s  for  this 
purpose  only. 
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Figure  3.  Hourly  solar  irradiance  values  at  5(fN,  145®W  on  Julian  date  10  using  the 
albedo  scheme  by  Coakley  (1979)  without  wind  speed  correction  -  solid;  and  with  wind 
speed  correction  -  dashed. 


C  SPECTRAL  DECOMPOSITION 

The  attenuation  with  depth  of  solar  irradiance  in  the  ocean  varies  considerably 
across  the  spectrum.  This  leads,  in  part,  to  the  blue  color  of  ocean  waters.  Thus  equation 
(1.1)  is  more  correctly  written 


13 


Hz)  = 


(2.10) 


where  k  is  wavelength.  To  approximate  this  equation,  the  solar  spectrum  can  be 
decomposed  into  a  finite  number  of  bands,  with  each  band  having  its  own  characteristic 
values  for  /o  and  k.  The  resulting  approximation  to  equation  (2.10)  is 

'fc)  =  ■  <2"> 

i»i 

Simpson  and  Dickey  (1981)  made  the  first  examination  of  the  effect  of  separate 
spectral  band  auenuation  versus  combined  attenuation  on  mixed  layer  dynamics.  In  their 
study,  the  spectral  decomposition  was  focused  on  the  infrared  spectrum.  The  visible  and 
ultraviolet  portion  of  the  solar  spectrum  were  divided  into  only  two  wavelength 
components,  200-600  nm  and  600-9(K)  nm.  Dynamic  studies  of  fresh  water  lakes  (e.g. 
Jassby  and  Powell,  1975),  have  also  used  spectral  considerations  in  the  representation  of 
solar  heat  flux. 

In  this  study,  decomposition  of  the  visible  spectrum,  as  well  as  the  ultraviolet  and 
infrared,  is  emphasized.  The  total  solar  spectrum  has  been  divided  into  1 1  bands  that 
encompass  the  ultraviolet,  visible,  and  infrared  portions  of  the  spectrum  as  described  in 
table  1 .  The  last  three  bands  constitute  the  infrared  region  of  the  solar  spectrum.  Values 
of  k  in  the  infrared  spectrum  are  very  large  relative  to  values  in  the  visible  bands  (A  - 
F).  indicating  the  infrared  portion  of  solar  irradiance  is  absorbed  very  near  to  the  ocean 
surface.  A  common  method  of  treating  this  region  of  the  spectrum  in  mixed  layer 
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m(xleling  is  to  assume  that  all  of  the  infrared  energy  is  absorbed  at  the  surface  The 
importance  of  more  accurately  including  the  slight  penetration  of  infrared  energy  in  a 
bulk  mixed  layer  m(xlel  is  examined  in  the  following  chapters.  The  results  of  Simpson 
and  Dickey  (1981)  indicate  that  consideration  of  this  penetration  'S  important. 


TABLE  1.  DEFINITION  OF  SPECTRAL 
BANDS. 


BAND 

WAVELENGTH(nm) 

UVl 

<  350 

UV2 

350  -  400 

A 

400  -  450 

B 

450  -  500 

C 

500  -  550 

D 

550  -  600 

E 

600  -  650 

F 

650  -  700 

IRl 

700  -  1000 

IR2 

1000  -  2000 

1R3 

>  2000 

The  /a  values  for  each  hand  of  the  spectrum  can  be  evaluated  through  use  of  a 


spectral  weighting  factor  a,  that  alters  equation  (2.11)  as  follows; 


(2.12) 


s 

nz)  « 

i  =  l 

The  a,  values  vary  with  F  anti  0.  Henderson  (1970)  indicates  that  slight  changes  in  the 
chemical  composition  of  the  atmosphere  significantly  alter  the  spectral  distribution  of 
sunlight,  especially  in  the  ultraviolet  and  infrared  regions.  Robinson  (1966)  analyzes  the 
relationship  between  atmospheric  turbidity,  which  represents  the  atmospheric  chemical 
content,  and  spectral  irradiance.  He  also  presents  a  table  that  relates  spectral  irradiance 
for  a  clean,  dry  atmosphere  to  values  of  (f>.  Robinson’s  table  was  used  in  this  study  to 
determine  the  values  of  a,  as  a  funcUon  of  <^.  Table  2  provides  values  for  each  15 
degrees  of  solar  altitude.  Inaccuracies  incurred  by  the  assumpdon  of  a  clean,  dry 
atmosphere  are  predominately  in  the  infrared  and  ultraviolet  bands  where  the  absorption 
effect  of  water  vapor  dominates.  Weighting  factors  for  discrete  values  are  linearly 
interpolated  from  tabulated  values.  Values  for  a,  in  the  visible  spectrum  were  compared 
with  surface  irradiance  measurements  of  Pak  et  al.  (1986)  conducted  as  part  of  the 
Optical  Dynamics  Experiment  (ODEX).  Errors  of  less  than  5  percent  were  achieved  in 
the  visible  bands  (A  -  F). 

Figure  4  provides  some  examples  of  the  variation  of  Ui  for  the  selected  bands.  The 
major  effect  of  a  decreasing  <(>  is  that  more  of  the  solar  energy  is  contained  in  the  longer 
wavelength,  or  "red"  bands,  as  more  energy  from  the  shorter  wavelengths  is  scattered  out 
of  the  illumination  path  by  Rayleigh  scattering  because  of  the  greater  path  length  at  low 
sun  angles.  This  leads  to  the  red  appearance  of  the  sun  at  low  altitudes. 
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TABLE  2.  WEIGHTING  FACTORS  {a,)  USED  TO  DETERMINE  THE  AMOUNT 
OF  TOTAL  INSOLATION  CONTAINED  IN  EACH  SPECTRAL  BAND  AS  A 
FUNCTION  OF  SOLAR  ALTITUDE  (<|>). 


SOLAR  ALTITUDE  (DEG) 

BAND 

0 

15 

30 

45 

60 

75 

90 

UVl 

.000 

.002 

.006 

.007 

.010 

.018 

.025 

UV2 

.000 

.012 

.018 

.024 

.029 

.035 

.043 

A 

.012 

.039 

.051 

.058 

.061 

.067 

.073 

B 

.021 

.058 

.067 

.071 

.074 

.077 

.080 

C 

.043 

.065 

.070 

.071 

.072 

.074 

.075 

D 

.050 

.069 

.071 

.071 

.071 

.071 

.070 

E 

.061 

.069 

.069 

.068 

.067 

.067 

.066 

F 

.064 

.067 

.065 

.064 

.063 

.062 

.060 

IRI 

.309 

.273 

.259 

.252 

.248 

.238 

.230 

IR2 

.343 

.275 

.255 

.248 

.242 

.229 

.220 

IR3 

.097 

.071 

.069 

.066 

.063 

.062 

.058 

The  sum  of  the  weighting  factors  in  bands  IRI,  1R2  and  IR3  represents  the  total 
percentage  of  insolation  that  is  in  the  infrared  portion  of  the  spectrum.  Ivanoff  (1977) 
and  Robinson  (1966)  present  additional  estimates  for  the  percentage  of  the  solar  spectrum 
in  the  infrared  region  for  clear  skies.  Figure  5  is  a  comparison  of  the  total  infrared 
irradiance  calculated  by  the  spectral  decomposition  used  in  this  study  with  those  two 
estima'tvjs.  It  should  be  noted  that  in  the  estimation  by  Robinson,  infrared  is  defined  as 
wavelengths  greater  than  740  nm,  while  Ivanoff  uses  700  nm  for  the  minimum  infrared 
wavelength  concordant  with  this  study.  This  explains  the  smaller  percentage  values 
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Figure  4.  Values  of  the  spectral  weighting  factor,  a,,  for  the  11  spectral  bands  for 
various  solar  altitudes.  Based  on  the  work  by  Robinson  (1966). 
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illustrated  in  Robinson’s  curve  at  all  <t>  values. 

It  should  be  noted  that  the  infrared  radiation  discussed  here  refers  only  to  that 
portion  of  the  solar  spectrum  and  not  to  the  infrared  radiation  emitted  by  atmospheric 
constituents.  Typically,  the  solar  spectrum,  including  the  infrared,  is  referred  to  (in 
oceanographic  literature)  as  "short-wave"  radiation,  while  atmospheric  radiation  is 
referred  to  as  "long- wave"  radiation.  This  terminology  is  derived  from  the  fact  that  the 
majority  of  atmospheric  emission  is  at  wavelengths  greater  than  7(X)0  nm,  or  what  is 
commonly  referred  to  as  the  "far  infrared. "  Solar  infrared  radiation  is  almost  exclusively 
in  the  "near  infrared"  (less  than  3000  nm).  Lind  and  Katsaros  (1982)  developed  a  model 
for  predicting  the  net  long- wave  radiation  for  various  atmospheric  conditions.  Long-wave 
radiation  is  not  considered  in  this  study. 

D.  REFRACTION  OF  DIRECT  SOLAR  INSOLATION 

The  total  downward  irradiance,  as  discussed  in  section  B  above,  can  be  divided 
into  direct  and  diffuse  components.  The  directionality  of  the  direct  comjionent  affects  not 
only  the  value  of  surface  reflectance  but  also  the  path  along  which  this  component 
penetrates  the  ocean.  If  the  sun  is  considered  as  a  point  source  of  solar  irradiance,  the 
direct  component  is  refracted  as  it  crosses  the  ocean  surface  (equation  2.7).  This 
refractive  effect  has  been  included  in  heat  flux  analyses  for  fresh  water  bodies  (e.g. 
Jassby  and  Powell,  1975),  but  not  in  similar  analyses  for  the  ocean.  The  impact  of  this 
effect  will  be  examined  in  the  following  chapters  based  on  the  representation  that  follows. 
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Figure  5.  Percent  of  the  total  insolation  contained  in  the  infrared  spectral  region  as 


calculated  by  Robinson  (1966)  -  solid;  Ivanoff  (1977)  -  dashed;  spectral  decomposition 


-  dotted. 


1 .  Determination  of  Direct  and  Diffuse  Insolation 

To  accurately  quantify  the  amount  of  radiation  being  refracted,  the  relative 
amount  of  direct  insolation  included  in  the  total  insolation  must  be  determined.  Jassby 
and  Powell  (1975)  performed  a  study  at  Castle  Lake  (Oregon)  during  a  period  when 
cloudless  skies  predominated.  They  therefore  assumed  that  all  insolation  was  direct. 
Iiradiance  measurements  in  the  region  indicated  that  this  assumption  was  valid  for  the 
high  altitude  Castle  Lake  region.  At  sea  level,  even  under  clear  skies,  the  diffuse 
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component  must  be  considered,  especially  in  the  shorter  wavelengths  due  to  the  great 
amount  of  atmospheric  scattering.  Curves  developed  by  Henderson  (1970)  based  on  work 
by  Hinzpeter  (1957)  prescribe  the  fractional  amount  of  diffuse  radiation  in  the  total 
insolation  as  a  function  of  wavelength  and  solar  altitude  for  a  Rayleigh  atmosphere 
(Figure  6).  A  Rayleigh  atmosphere  includes  only  scattering  by  molecular  components  the 
size  of  which  are  much  less  than  the  wavelength  of  the  electromagnetic  energy  being 
scattered;  scattering  by  particulates  is  not  considered.  The  diffuse  component  of  all 
wavelengths  increases  as  <t>  decreases,  and  the  diffuse  component  is  relatively  large  for 
shorter  wavelengths  at  all  solar  altitudes.  The  curves  in  Figure  6  were  used  to  construct 
the  values  for  a  fractional  coefficient  (B)  that  prescribes  the  amount  of  total  irradiance 
contained  in  the  direct  and  diffuse  components  as  follows: 


^0  direct  ~  P  4 

(2.13) 

^diffuse  (1  ■  P)A) 

(2.14) 

Table  3  lists  B  values  for  each  ten  degrees  of  <^,  similar  to  the  weighting  factors  listed  in 
Table  2. 

2.  Refractive  Effects  of  the  Ocean 

The  actual  path  over  which  the  absorption  of  direct  solar  irradiance  occurs  in 
the  ocean  is  prescribed  as  a  function  of  6,  the  angle  of  refraction  measured  from  the 
horizontal  as  follows: 
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Wortte<rg(^  (pmi 


Figure  6.  Curves  produce  by  Henderson  (1970)  illustrating  the  percentage  of  total 
insolation  contained  in  the  diffuse  (Sky  radiation)  and  direct  (Solar  radiation) 
components  as  a  function  of  solar  altitude  and  wavelength. 


TABLE  3.  FRACTIONAL  COEFFICIENT  (B)  VALUES  USED  TO  DETERMINE 
THE  AMOUNT  OF  IRRADIANCE  CONTAINED  IN  THE  DIRECT  AND  DIFFUSE 
COMPONENTS  OF  INSOLATION  FOR  EACH  SPECTRAL  BAND  AS  A 
FUNCTION  OF  SOLAR  ALTITUDE. 


SOLAR  ALTITUDE  (DEG) 

BAND 

0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

UVl 

1.0 

.98 

.87 

.80 

.69 

.63 

.60 

.57 

.55 

.53 

UV2 

1.0 

.94 

.74 

.60 

.52 

.46 

.43 

.40 

oo 

.37 

A 

1.0 

.83 

.60 

.45 

.39 

.35 

.32 

.30 

.28 

.28 

B 

.99 

.73 

.48 

.35 

.28 

.25 

.23 

.21 

.21 

.20 

C 

.95 

.60 

.40 

.28 

.23 

.20 

.18 

.17 

.16 

.16 

D 

.90 

.51 

.33 

.23 

.18 

.16 

.14 

.14 

.13 

.13 

E 

.78 

.42 

.26 

.19 

.16 

.14 

.12 

.12 

.11 

.11 

F 

.65 

.38 

.23 

.16 

.13 

.11 

.10 

.09 

.09 

.09 

IRl 

.50 

.22 

.12 

.09 

.08 

.07 

.07 

.06 

.06 

.06 

IR2 

.41 

.08 

.05 

.04 

.04 

.03 

.03 

.02 

.02 

.02 

IRS 

.24 

.07 

.05 

.03 

.03 

.02 

.02 

.02 

.02 

.02 

sin  6 


(2.15) 


Using  Snell’s  law,  equation  (2. 15)  can  be  rewritten  in  terms  of  <!>  and  the  relative  index 
of  refraction  (n): 


s  = 


(2.16) 


\Jl  -  cos^6  y/l  -  /i'^cos^4) 

The  total  irradiance  beneath  the  surface  cannot  be  fully  represented  by  separate 
direct  and  diffuse  components  (Jerlov  .1976).  As  depth  increa.ses.  some  of  the  direct 
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component  is  scattered  by  seawater  constituents  and  becomes  diffuse.  Thus  the  direct 
component  is  attenuated  at  a  rate  greater  than  that  predicted  by  exponential  attenuation, 
and  the  diffuse  component  is  attenuated  less.  The  magnitude  of  this  effect  is  a  function 
of  the  seawater  turbidity  with  additional  spectral  implications.  Surface  waves  cause  a 
departure  from  Snell’s  law'  by  inducing  fluctuations  in  the  direction  of  the  refracted  direct 
component.  This  effect  will  be  neglected  in  this  study,  and  two-component  penetration 
with  the  use  of  Snell’s  law  provides  an  accurate  representation  of  the  near-surface 
irradiance.  Since  this  is  the  region  where  the  majority  of  the  downward  irradiance  is 
absorbed,  this  method  should  not  induce  noticeable  error  into  the  analyzed  effects  on 
mixing  dynamics  and  thermal  stratification. 

Combining  equations  (2-13))  through  (2-16),  the  total  downward  irradiance 
for  a  specific  spectral  band  l,{z)  can  be  calculated  by  the  following  combination  of  direct 
and  diffuse  downward  irradiance: 


k,z 


(2.17) 


Combining  the  parameterizations  of  spectral  decomposition  and  refraction, 
equation  (1.1)  can  now  be  written  as; 
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(2.18) 


m 


11 


/oE 


P.)e 


yfl  -  n  “oos?^  j 


*,z 


where  the  sum  is  over  the  1 1  spectral  bands  described  in  section  C. 

If  equation  (1.1),  the  simple  exponential  with  a  single  k,  is  used  to  describe 
the  penetration  of  sunlight  into  the  ocean,  a  fraction  of  U  representing  the  infrared  solar 
radiation  must  be  considered  as  being  absorbed  exactly  at  the  ocean  surface.  The  value 
fraction  commonly  used  in  mixed  layer  models  is  1/2.  Equation  (1.1)  therefore  becomes 

/(c)  -  (2.19) 

for  actual  use  in  mixed  layer  models  such  as  Garwood  (1977)  and  others. 

E.  ATTENUATION  COEFFICIENTS 

Evaluation  of  the  above  equations  requires  that  accurate  attenuation  coefficient 
values  be  obtained  for  the  specific  cases  examined.  The  single  k  used  in  equation  (1.1) 
describes  only  the  visible  portion  of  the  solar  spectrum  and  is  usually  estimated  from 
Secchi-disk  observations.  The  infrared  portion  of  the  spectrum  is  assumed  to  be  absorfied 
exactly  at  the  surface  and  attenuation  of  the  ultraviolet  portion  is  commonly  neglected. 
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A  commonly  used  algorithm  for  obtaining  k  from  Secchi-disk  observations  is  /t  =  \.lld, 
where  d  is  the  Secchi  depth  in  meters  (Myers  et  al.,  1969).^ 

Direct  measurement  of  irradiance  attenuation  with  a  pyranometer  provides  the  most 
accurate  source  of  spectral  k,  values.  This  type  of  irradiance  measurement  allows 
calculation  of  k,  values  as  a  function  of  depth  as  well  as  wavelength.  The  value  of  k 
varies  with  depth  due  to  layering  of  absorptive  particulates  such  as  phytoplankton,  and 
the  nonlinear  effect  of  surface  waves  on  ocean  optical  properties.  Vertical  variation  of 
k,  should  be  considered  in  detailed  studies  of  ocean  optical  conditions  at  specific  locations 
where  measurements  are  made. 

Calculation  of  k,  value';  from  colorimetric  satellite  data  is  another  method.  Au.stin 
(1981)  prescribes  algv>i  -  .>ms  for  calculating  k,  for  different  wavelengths  using  data  from 
the  Coastal  Zone  Color  Scanner  (CZCS).  Future  colorimeters  will  incorporate  more 
spectral  ^'hannels  and  will  provide  an  accurate  and  readily  available  source  of  k, 
measurements.  The  k,  values  from  satellite  observations  represent  a  vertical  mean  value, 
as  variation  with  depth  can  not  be  measured. 

Jerlov  (1976)  constructed  an  ocean  atlas  of  optical  water  types  and  provides  k 
values  for  each  water  type  in  the  ultraviolet  and  visible  spectrum  based  on  many  years 
of  spectral  irradiance  observations.  The.se  optical  water  types  provide  a  baseline  of  values 

’The  value  1 .7  is  due  to  Poole  and  Atkins  (1929),  who  based  this  value  on  only  14 
different  data  points  -  their  work  has  been  accepted  uncritically  by  biologists  (mainly). 
It  is  more  accurate  to  consider  k  =  xld  where  x  is  a  value  ranging  from  1.0  to  2.0, 
depending  on  the  water.  (Walker,  1980,  recalculated  the  data  of  Poole  and  Atkins  and 
found  X  =  1 .45  to  be  more  accurate.)  See  Preisendorfer  (1986)  regarding  caveats  on  the 
use  of  Secchi  data. 
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when  direct  measurements,  performed  either  in  situ  or  remotely,  are  unavailable.  Remote 
measurements  of  a  few  spectral  k,  values  may  also  identify  an  oceanic  region  in  question 
as  a  specific  Jerlov  water  type,  allowing  use  of  the  Jerlov  ki  values  at  wavelengths  that 
are  not  measured  by  the  colorimeter.  The  k,  values  in  the  infrared  spectrum  are 
extremely  high,  with  all  infrared  radiation  being  absorbed  in  the  first  few  centimeters. 
Measurements  indicate  that  particulates  in  seawater  do  not  play  an  important  role  as  the 
pure  water  molecules  are  responsible  for  the  majority  of  absorption  iii  dit  infrared. 
Defant  (1961)  provides  typical  values  for  ki  at  infrared  wavelengths  that  are  used  in  this 
study  for  all  water  types.  Table  4  lists  k,  values  for  the  five  open  ocean  Jerlov  water 
types,  and  Figure  7  provides  vertical  profiles  of  irradiance  (/)  for  the  same  water  types. 
Equation  (2.18)  is  used  to  calculate  the  irradiance  profiles  in  Figure  7,  with  a  value  of 
5(X)  W/m’  for  /o;  solar  altitude  is  assumed  to  be  90  degrees. 

The  k,  values  in  the  vicinity  of  (Kean  station  "P"  are  of  particular  interest  in  this 
study,  as  numerical  experiments  in  chapter  IV  are  performed  using  meteorological  and 
(Keanographic  observations  performed  by  vessels  at  (Kean  station  "P"  from  1955  through 
1969.  Jerlov  describes  the  (Kean  waters  in  the  vicinity  of  "P"  as  water  type  11.  but 
irradiance  measurements  (Austin,  1991)  provide  k,  values  that  correspond  to  type  IB 
waters.  Values  of  k,  used  in  the  numerical  evaluation  of  chapter  IV  represent  type  IB 
waters.  Selection  of  a  representative  k  for  the  single  exponential  scheme  of  equation 
(1.1)  is  discussed  in  chapter  111. 

Optical  properties  of  (Kean  waters  vary  with  changes  in  the  amount  of  dissolved 
solids  and  particulates  in  the  water.  Phytoplankton  bkxims  and  river  run-off  cause  such 
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TABLE  4.  SPECTRAL  ATTENUATION  COEFFICIENT  k,  VALUES  (m ')  FOR 
THE  FIVE  JERLOV  OPEN  OCEAN  WATER  TYPES  (JERLOV,  1976;  DEFANT, 
1961) 


JERLOV  WATER  TYPE 

BAND 

I 

lA 

IB 

11 

III 

UVl 

.106 

.129 

.160 

.272 

.485 

UV2 

.038 

.052 

.066 

.122 

.220 

A 

.022 

.031 

.042 

.081 

.160 

B 

.018 

.025 

.033 

.062 

.116 

C 

.043 

.048 

.054 

.076 

.116 

D 

.089 

.094 

.099 

.115 

.148 

E 

.305 

.310 

.315 

.335 

.375 

F 

.420 

.430 

.435 

.465 

.520 

IRl 

10.242 

10.242 

10.242 

10.242 

10.242 

IR2 

65.669 

65.669 

65.669 

65.669 

65.669 

IR3 

4620.2 

4620.2 

4620.2 

4620.2 

4620.2 

changes.  Variation  in  the  amount  of  phytoplankton  is  the  primary  mechanism  responsible 
for  optical  variation  in  colder,  open  ocean  waters  such  as  those  surrounding  ocean  station 
"P."  These  variations  occur  periodically  as  seasonal  changes  in  water  temperature, 
mixing,  and  sunlight  availability  influence  phytoplankton  production.  Secchi  observations 
over  a  period  of  years  can  be  examined  to  reveal  seasonal  fluctuations  in  optical 
characteristics.  Figure  8  illustrates  Secchi  measurements  in  the  vicinity  of  ocean  station 
"P"  based  on  the  date  of  measurement.  The  data  set  used  for  Figure  8  was  provided  by 
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the  National  Oceanographic  Data  Center  (NODC),  based  on  recorded  oceanographic 
measurements  since  1960.  A  strong  seasonal  periodicity  in  Secchi  measurements  cannot 
be  inferred  from  Figure  8;  however,  a  trend  of  clearer  water  during  the  winter  and  spring 
months  is  evident  as  phytoplankton  production  is  strongest  in  late  summer  and  fall  after 
the  period  when  the  daily  amount  of  insolation  is  maximum. 


Figure  7.  Vertical  uradiance  profiles  for  the  5  Jerlov  open  ocean  w-ater  types. 
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Figure  8.  Secchi  observations  in  the  vicinity  of  ocean  station  "P"  grouped  by  year 
day  of  measurement. 


in.  ANALYTICAL  EVALUATION 


The  algorithms  developed  in  Chapter  II  are  applied  to  expressions  that  relate  the 
solar  heat  ttux  in  the  upper  ocean  to  the  dynamic  processes  responsible  for  its  density 
structure.  The  production  and  dissipation  of  turbulent  kinetic  energy  and  the  associated 
thermodynamic  processes  include  the  penetration  of  solar  radiation  in  the  determination 
of  both  depth  and  temperature  of  the  ocean  mixed  layer.  The  sensitivity  of  these 
prtKesses  to  the  various  parameterization  schemes  discussed  in  chapter  II  is  examined. 
The  three  schemes  are  defined  as  follows  for  future  reference: 

Scheme  1  -  Simple  exponential  attenuation  of  irradiance  using  a 

single  constant  k  value  lequation  (2.19)]. 

Scheme  2  -  Exponential  attenuation  of  irradiance  in  individual 

spectral  bands  using  specific  A  values  for  each  band  that 
vary  with  .solar  altitude  [equation  (2.121). 

Scheme  ,1  -  Exponential  attenuation  of  individual  spectral  bands 

combined  with  refractive  effects  on  direct  sunlight 
[equation  (2.18)]. 

Regarding  notation  of  physical  quantities  such  as  temperature,  density,  and  velocity 
components,  an  over-bar  denotes  the  en.semble  mean  of  a  quantity  and  a  prime  denotes 
the  fluctuating  part,  for  example 

T  =  T  T'  and  w  =  w  +  w' 
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A.  TURBULENT  ENTRAINMENT 


1 .  Effective  Surface  Buoyancy  Flux 

Following  the  turbulence  closure  scheme  of  Garwood  (1977)  which 
incorporates  bulk  buoyancy  and  momentum  equations,  the  budget  for  the  vertical 
component  of  turbulent  kinetic  energy  can  be  written  as 


»«,((£)  -3{w'w'»(£)'«- 


m. 


~fh 


m, 


{E) 


(3.1) 


where  h  is  mixed  layer  depth;  £  =  (u^  +  +  w^)  is  proportional  to  the  turbulent  kinetic 

energy;  /w,,  m.,  and  Wj  are  dimensionless  constants  of  proportionality;  uJ),  is  the 
effective  downward  surface  buoyancy  flux;  h  is  buoyancy  defined  as: 


where  p  is  defined  by 


gipQ-p) 

Po 


(3.2) 


p  =  Po[  1  -  a(e  -  Oq)  +  P(s  -  So)] .  (2.3) 

In  (3.2)  and  (3.3),  g  is  gravity,  p  density,  s  salinity,  T  temperature,  a  and  B  the 
expansion  coefficients  for  heat  and  salt,  respectively.  The  subscript  zero  denotes 
representative  reference  values.  In  this  study,  only  the  relationship  between  temperature 
and  buoyancy  are  of  interest,  thus  s  =  So,  and  changes  in  density  are  caused  only  by 
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III.  ANALYTICAL  EVALUATION 


The  algorithms  developed  in  Chapter  11  are  applied  to  expressions  that  relate  the 
solar  heat  flux  in  the  upper  ocean  to  the  dynamic  processes  responsible  for  its  density 
structure.  The  production  and  dissipation  of  turbulent  kinetic  energy  and  the  associated 
thermodynamic  processes  include  the  penetration  of  solar  radiation  in  the  determination 
of  both  depth  and  temperature  i)f  the  ocean  mixed  layer.  The  sensitivity  of  these 
pn.)cesses  to  the  various  parameterization  schemes  discussed  in  chapter  II  is  examined. 
The  three  schemes  are  defined  as  follows  for  future  reference: 

Scheme  1  -  Simple  exponential  attenuation  of  irradiance  using  a 

single  constant  k  value  [equation  (2.19)]. 

Scheme  2  -  Exponential  attenuation  of  irradiance  in  individual 

spectral  bands  using  specific  k,  values  for  each  band  that 
vary  with  .solar  altitude  [equation  (2.12)[. 

Scheme  3  -  Exponential  attenuation  of  individual  spectral  bands 

combined  with  refractive  effects  on  direct  sunlight 
[equation  (2. 18)[. 

Regarding  notation  of  physical  quantities  such  as  temperature,  density,  and  velocity 
comptments,  an  over-bar  denotes  the  en.semble  mean  of  a  quantity  and  a  prime  denotes 
the  fluctuating  part,  for  example 

T  T  +  T'  and  w  =  w  +  w' 
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A.  TURBULENT  ENTRAINMENT 


1 .  Effective  Surface  Buoyancy  Flux 

Following  the  turbulence  closure  scheme  of  Garwwxl  (1977)  which 
incorporates  bulk  buoyancy  and  momentum  equations,  the  budget  for  the  vertical 
component  of  turbulent  kinetic  energy  can  be  written  as 


\t  /  + 


m. 


m 


1  ) 


(E) 


(3.1) 


where  h  is  mixed  layer  depth;  E  =  (u^  +  +  h^)  is  proportional  to  the  turbulent  kinetic 

energy;  m.,  and  w,  are  dimensionless  constants  of  proportionality;  uJ?.  is  the 
effective  downward  surface  buoyancy  flux;  h  is  buoyancy  defined  as: 


where  p  is  defined  by 


^(pQ-  P) 
Po 


(3.2) 


p  -  Po[  1  -  a(0  -  %)  +  ^s  -  So)]  •  (2-3) 

In  (3.2)  and  (3.3),  g  is  gravity,  p  density,  s  .salinity,  T  temperature,  a  and  B  the 
expansion  coefficients  for  heat  and  salt,  re.spectively.  The  sub.script  zero  denotes 
representative  reference  values.  In  this  study,  only  the  relationship  between  temperature 
and  buoyancy  are  of  interest,  thus  s  =  Sq,  and  changes  in  density  are  caused  only  by 
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temperature  fluctuations.  An  approximation  of  local  horizontal  homogeneity  with  regard 
to  mean  turbulence  fields  is  used,  allowing  horizontal  turbulence  advection  to  be 
neglected.  Vertical  advection  of  turbulence  via  the  vertical  component  of  turbulent 
velocity  is  therefore  the  principal  mechanism  causing  entrainment.  This  one-dimensional 
analysis  allows  the  role  of  optical  effects  to  be  emphasized. 

The  effective  downward  surface  buoyancy  flux,  uJ?.,  is  used  to  prescribe  the 
effect  of  surface  heat  flux  and  solar  radiation  on  the  amount  of  turbulent  kinetic  energy 
in  the  fully  turbulent  bulk  mixed  layer  and  is  defined  by  the  following  equation; 

r°(rV)*  +  (3.4) 

The  value  of  u.h.  may  be  computed  using  thermodynamic  relationships  in  the  mixed 
layer.  The  turbulent  heat  flux  (Tw')  can  be  analyzed  using  the  heat  budget  equation  for 
the  mixed  layer: 

^  ,  _Q_  (3.5) 

dt  dz  pCp 

where  Q  is  the  vertical  divergence  of  irradiance  (dlldz),  which  represents  absorption  of 
solar  radiation,  as  described  in  chapter  1,  and  Cp  is  the  specific  heat  of  seawater  at 
constant  pressure,  approximately  equal  to  3990  J  kg  '  K  ‘  (Gill,  1982). 

The  mixed  layer  temperature  is  assumed  constant  with  respect  to  z,  thus  the 


left  and  right  hand  sides  of  equation  (3.5)  are  constant  in  z  and  can  be  rewritten  as; 


il  -  A 

dt  "  dz 


-{t'w') 


constant  in  z. 


(3.6) 


The  sum  inside  the  brackets  of  equation  (3.6)  is  therefore  linear  with  respect 
to  z  and  can  be  expressed  as; 


T'w 


pc. 


(rv'),  ^ 


pc. 


(■  i)  ■  r ■  7^ 


(3.7) 


Rearranging  equation  (3.7)  and  substituting  into  equation  (3.1),  uJ).  can  now 
be  defined  by  the  following  equation; 


ub 


^  r  J-di 

h  pc 


_  lag  ra 
h  J  -h 


{T'wX  - 
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1 


-  dt 
h' 


+  20^/*° 
h  J  -h 


which  can  be  simplified  to; 


pc. 


ag 


{rw')_ 


(3.8) 


^  -  i/W 


pc^h  J  -A  '  h''  i 


dz 


(3.9) 


where  (b  wOo  is  the  net  buoyancy  flux  at  the  surface,  due  to  the  net  heat  flux,  and  can  be 
represented  by; 


(b'w'\  =  ag(T'w')^  (3.10) 

The  term  {Tw)„  represents  the  net  sum  of  the  evapt>rative.  sensible,  latent  and  long  wave 
radiative  heat  fluxes  at  the  surface. 
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A  positive  value  for  uJ).  indicates  a  downward  buoyancy  flux  into  the  mixed 
layer,  caused  only  by  negative  heat  fluxes  since  salinity  fluxes  are  being  neglected.  The 
net  surface  buoyancy  flux,  (positive  upward)  affects  the  amount  of  turbulence  in 
the  mixed  layer.  If  net  heat  is  added  to  the  surface  less  than  zero],  a  statically 

stabilizing  effect  results  in  the  mixed  layer  due  to  induced  stratification  by  the  warming 
of  surface  waters.  Turbulent  energy  is  expended  to  counteract  this  stratification  process, 
reducing  the  amount  of  turbulent  energy  available  for  entrainment.  If  net  heat  is  removed 
from  the  surface,  static  in.stability  enhances  the  amount  of  turbulent  energy  available  for 
entrainment.  The  absorption  of  penetrating  solar  radiation  also  affects  the  amount  of 
turbulence  available  for  entrainment.  The  buoyancy  flux  at  the  surface  is  reduced  by  the 
solar  radiation  penetrating  below  the  surface,  which  represents  a  heat  or  buoyancy  flux 
into  the  water  below  the  surface.  If  absorption  occurs  predominately  near  the  top  of  the 
mixed  layer,  a  stabilizing  effect  a^ain  reduces  the  effective  turbulence  in  the  mixed  layer. 
However,  if  absorption  of  radiation  is  constant  throughout  the  mixed  layer,  the  effective 
turbulence  is  unaffected  by  solar  absorption. 

The  ef  fect  of  solar  radiation  on  turbulent  kinetic  energy  through  uJ).  can  be 
resolved  from  equation  (3.9)  by  defining  a  function  F  that  depicts  the  absorption  curve 
in  the  mixed  layer  as  follows: 

u.b,  =  -(^%  -  -R)F  .  (3.11) 
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2.  Application  of  Optical  Parameterizations 

The  functional  form  of  F  is  expressed  differently  depending  upon  the  scheme 
used  to  prescribe  the  attenuation  of  I.  The  following  three  equations  specify  F  for  each 
of  the  three  optical  schemes,  where  the  subscript  denotes  the  corresponding  scheme: 


1  ,  I[i  -  A(i 

2  2[  kh 


(3.12) 
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These  three  variations  for  Fean  be  compared  by  allowing  h  to  vary.  Figure 
9  represents  such  a  comparison.  The  k,  values  used  to  construct  Figure  9  were  selected 
for  Jerlov  type  IB  waters  as  described  in  chapter  II.  These  values  correspond  to  the 
waters  in  the  vicinity  of  ocean  statiiin  "P." 
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A  greater  than  zero  value  for  the  function  F  will  always  act  to  reduce  the 
predicted  amount  of  wind  generated  turbulence  in  the  mixed  layer,  as  solar  radiation  is 
absorbed  predominately  near  the  surface  and  thus  tends  to  stabilize  the  mixed  layer 
through  thermal  stratification,  while  the  turbulent  mixing  counteracts  this  stratification 
process.  This  is  illustrated  by  the  range  of  F  in  Figure  9.  For  a  deep  mixed  layer  (large 
h),  all  F  values  converge  toward  1 .0  as  nearly  all  solar  radiation  is  absorbed  at  relatively 
shallow  depths,  and  the  stabilizing  effect  is  maximum.  For  a  very  shallow  mixed  layer, 
a  significant  amount  of  absorption  may  occur  below  the  mixed  layer,  and  the  stabilizing 
effect  is  minimal.  Figure  9  illustrates  that  the  value  of  F  changes  considerably  for  the 
different  schemes  used  in  this  study,  thus  having  a  significant  impact  on  uJ).  and  the 
amount  of  turbulent  kinetic  energy  available  for  entrainment  at  the  base  of  the  mixed 
layer. 


B.  EFFECT  ON  TEMPERATURE  DISTRIBUTION 

The  temperature  of  the  mixed  layer  evolves  as  a  result  of  surface  heat  fluxes, 
entrainment,  and  horizontal  thermal  advection.  Advective  effects  are  neglected  in  this 
one-dimensional  analysis.  The  evolution  of  temperature  in  the  mixed  layer  can  therefore 
be  prescribed  by  integrating  equation  (3.5)  over  the  mixed  layer  and  recalling  that  dTIdt 
is  assumed  constant  with  respect  to  z  in  the  mixed  layer.  The  resulting  equation  is 


dT. 


ml 


dt  h 

L 

where  T„,  is  the  mixed  layer  temperature 
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Figure  9.  Values  of  F,  -  solid.  Fj  -  dashed,  and  Fj  -  dotted,  for  Jerlov  type  IB  water 


as  a  function  of  mixed  layer  depth  (h)  at  various  solar  altitudes. 


A  simple  experiment  is  performed  to  examine  the  optical  effect  on  mixed  layer 
temperature  change  It  is  assumed  that  the  net  surface  heat  flux  is  in  perfect 
thermodynamic  balance  KFhOo  =  0|,  the  entrainment  heat  flux  is  zero,  and  there  is  no 
change  in  mixed  layer  depth.  If  these  conditions  are  maintained  over  a  period  of  time, 
the  absorption  of  solar  radiation  provides  the  only  heat  input  to  the  mixed  layer,  and 
equation  (3.15)  becomes 


38 


(3.16) 


Values  for  /„  and  4>  are  calculated  for  Julian  dates  191-200  at  a  latitude  of  50°N.  The  k 
value  for  scheme  1  produces  an  equal  value  of  /h  for  all  three  schemes  at  a  solar  altitude 
of  90“.  (For  4>  ^  90“  there  is  no  refraction.)  This  would  make  the  solar  absorption  and 
corresponding  temperature  increase  equivalent  for  all  three  schemes  if  4>  was  held 
constant  at  90“.  Figure  10  therefore  accurately  reflects  the  sensitivity  of  mixed  layer 
temperature  evolution  to  the  variation  of  schemes  2  and  3  with  </>.  This  variability  with 
<i)  contributes  to  the  improved  accuracy  of  these  schemes  over  scheme  1 .  The  value  k  = 
0.08864  is  also  used  for  scheme  1  in  the  numerical  evaluation  of  chapter  IV. 

Using  equation  (3.16),  a  function  G  can  be  used  to  specify  the  relative  amount  of 
the  surface  irradiance  that  is  abst)rbed  in  the  mixed  layer  and  therefore  responsible  for 
the  V  arming  of  the  mixed  layer  in  a  simple  mrxlel: 
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Figure  10.  Increase  in  mixed  layer  temperature  as  a  result  of  different  solar 
absorption  parameterizations.  Solid  -  scheme  1;  dashed  -  scheme  2;  dotted  -  scheme 
3. 


^  =  Js-G  (3.17) 

dt  hpc^ 

Three  different  G  functions,  similar  to  the  F  functions  used  above  to  evaluate  the 
effective  surface  buoyancy  flux,  represent  each  of  the  attenuation  schemes.  The  different 
G  functions  are  defined  as  follows: 
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Gj  =  1  -O.Se*'’  (3.18) 

G.  =  1  -  (3-19) 

1=1 

-k,h 

11  — 

^3  =  1- 

1  =  1 

The  value  of  G  represents  the  fractional  amount  of  the  total  surface  irradiance  that  is 
absorbed  in  the  mixed  layer.  Values  of  G ,  G  ,  and  G,  as  a  function  of  4>  for  Jerlov 
water  type  IB  with  an  h  of  10  m  are  illustrated  in  Figure  12.  As  expected,  more  of  the 
surface  irradiance  is  absorbed  at  low  solar  altitudes  for  Gj  due  to  the  increased 
concentration  of  the  incident  solar  spectrum  into  the  longer  wavelengths  with  higher  k, 
values.  The  increased  absorption  due  to  refraction  is  demonstrated  by  the  difference 
between  G^  and  Gx. 

In  consideration  of  actual  mixed  layer  processes,  mixing  and  temperature  evolution 
can  not  be  separated,  as  the  two  roles  are  interdependent.  The  results  of  the  simple 
analyses  above  could  lead  to  the  hypothesis  that,  over  a  period  of  time,  the  mean  solar 
attenuation  can  be  accurately  represented  by  the  single  exponential  of  scheme  1  if  an 
accurate  k  is  used.  However,  the  accuracy  of  the  attenuation  scheme  during  short  time 
periods  may  be  crucial  to  accurate  mixed  layer  predictions,  as  individual  mixed  layer 
deepening  and  shallowing  events  are  triggered  or  enhanced  by  the  interaction  of 
atmospheric  forcing  and  heat  fluxes  over  small  time  scales.  The  overall  effect  of  these 


(3.20) 
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Values  of  "G"  Functions  (uniticss) 


different  schemes  on  mixed  layer  processes  is  examined  in  the  numerical  evaluation 
performed  in  chapter  IV. 


Figure  11.  Values  of  G,,  G2,  and  G3  as  a  function  of  <t>  for  Jerlov  water  type  IB. 
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IV.  NUMERICAL  EVALUATION 


In  this  chapter,  a  mixed  layer  model  is  used  to  examine  the  effects  of  the  optical 
parameterizations  developed  in  chapter  11.  Model  predictions  describing  the  upper  ocean 
thermal  structure  are  compared  to  examine  the  sensitivity  to  the  different  optical  schemes. 

A.  DESCRIPTION  OF  MODEL  AND  INITIALIZATION  DATA 

The  model  used  for  this  study  is  the  Garwood  (1977)  one-dimensional,  second  order 
closure,  bulk,  mixed  layer  model.  Horizontal  advection  is  not  included  in  the  model. 
Model  predictions  of  ocean  temperature  are  produced  by  the  model  for  the  upper  200 
meters  of  the  ocean  u.sing  a  one-meter  vertical  grid  resolution  and  first  order,  time- 
forward  finite  differencing.  Predictions  are  carried  out  hourly. 

Atmospheric  forcing  is  derived  from  meteorological  observations  performed  at 
ocean  station  "P"  (.'iO^’N,  145'’W)  during  the  years  1955-1969.  Values  for  total  insolation 
are  derived  for  t)cean  station  "P"  using  the  method  of  Reed  (1977),  as  described  in 
chapter  11.  Long  wave  radiation  is  calculated  from  the  method  of  Reed  and  Halper 
n  (1975).  Unless  otherwise  mentioned,  cloudless  conditions  are  used  in  the  calculation 
of  insolation  and  long  wave  radiation.  The  assumption  of  cloudless  skies  for  continuous 
periods  of  time  is  not  repre.sentative  of  actual  conditions  in  the  vicinity  of  ocean  station 
"P,"  but  this  assumption  en.sures  accuracy  in  the  u.se  of  a,  and  B,  values  for  schemes  2  and 
.1.  Cloudless  skies  also  provide  a  most-sensitive  scenario,  as  the  variation  of  schemes  2 
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and  3  from  scheme  1  is  maximum  for  clear  skies  when  the  direct  sunlight  component, 
which  varies  with  4),  is  maximum.  The  intent  of  this  evaluation  is  to  compare  the  results 
of  the  different  schemes,  not  to  verify  which  scheme  provides  the  most  accurate 
predictions.  To  accurately  verily  the  model  using  the  different  schemes,  realistic 
atmospheric  effects  (cloud  cover,  water  vapor  and  aerosols)  must  be  incorporated.  Values 
of  the  empirical  constants  used  for  physical  processes  not  covered  in  this  study  will  also 
have  to  be  re-evaluated. 

The  initial  vertical  temperature  structure  for  each  prediction  period  is  derived  from 
measurements  at  ocean  station  "P"  during  1955-1966.  As  described  in  chapter  11,  kj 
values  in  the  ultraviolet  and  visible  bands  represent  Jerlov  (1976)  water  type  IB,  while 
values  in  the  infrared  are  assumed  constant  for  all  water  types. 

B.  COMPARISON  OF  ALBEDO  SCHEMES 

The  mtxlel  is  used  to  compare  the  use  of  a  constant  daily  mean  albedo  with  a  value 
that  varies  hourly  as  a  function  of  the  mean  solar  altitude  (Coakley,  1979)  and  wind  speed 
(KaLsaros  et  al. ,  1985)  during  the  hour.  A  typical  lO-day  period  during  summer,  when 
lo  is  maximum,  is  used  to  test  the  albedo  effect.  Figures  12A  and  12B  provide  results 
v)f  this  comparison,  where  constant  albedo  is  represented  by  the  sr)lid  curves  and  the 
varying  albedo  by  dashed  curves. 

Figure  12A  indicates  that  the  varying  albedo  scheme  provides  more  heat  from  solar 
radiadon  into  the  ocean  during  the  day.  This  result  is  expected,  as  the  varying  albedo 
is  less  than  the  daily  mean  albedo  when  the  sun  is  at  high  altitudes,  where  h  is  large,  and 
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greater  than  the  mean  at  low  altitudes,  when  lo  is  less.  The  value  of  uJ).  in  Figure  12A 
indicates  increased  solar  absorption  near  the  surface  for  the  varying  albedo  during  periods 
of  high  (l>.  Consequently,  the  amount  of  turbulent  entrainment  during  deepening  events 
is  reduced.  This  results  in  a  shallower  mixed  layer,  also  reflected  in  Figure  12A. 
Greater  temperature  increase  is  depicted  in  Figure  12B  throughout  the  period  due  to 
increased  solar  absorption  at  all  depths  for  the  varying  albedo. 

C.  THE  ROLE  OF  ABSORPTION 

Before  applying  the  three  optical  parameterization  schemes  to  the  mixed  layer 
model,  vertical  profiles  of  solar  attenuation  and  absorption  are  examined.  For  schemes 
2  and  3,  the  values  of  a,  and  6,  in  equations  (2.12)  and  (2.18)  vary  as  a  function  of  <^, 
resulting  in  a  variation  of  the  attenuation  and  absorption  with  <^.  Figures  13A  and  13B 
provide  attenuation  and  absorption  profiles,  respectively,  for  different  values  of 

Figure  13A  illustrates  the  greater  attenuation  of  /  at  shallow  depths  for  schemes  2 
and  3,  especially  at  lower  4>  values  when  more  solar  radiation  is  concentrated  in  the 
longer  wavelength  bands  with  higher  k,  values.  The  increased  role  of  refraction  for  low 
altitudes  is  apparent  in  the  attenuation  profiles  for  scheme  3.  The  profile  of  Q  in  Figure 
1 3B  represents  the  relative  amount  of  solar  radiation  transferred  to  the  ocean  waters  as 
heat.  It  is  apparent  that  for  schemes  2  and  3,  more  warming  should  occur  at  the  very 
shallow  and  deep  regions,  while  less  warming  occurs  at  the  intermediate  depths.  The 
effect  of  these  inflection  points  in  the  absorption  profiles  are  important  factors  in  the 
response  of  the  mrxlel  to  the  individual  schemes. 
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Figure  12A.  Atmospheric  forcing  and  mixed  layer  depth  time  series  for  constant 


(solid)  and  varying  (dashed)  albedo  during  summer  conditions. 


46 


Temperature  ((leg  C)  Temperature  (dep  C)  Temperature  (deg  C)  Temperature  (deg  C) 


Figure  12B.  Vertical  temperature  structure  time  series  for  constant  (solid)  and 


varying  (dashed)  albedo  during  summer  conditions. 
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Figure  13A.  Vertical  attenuation  profiles  representative  of  waters  in  the  vicinity  of 
ocean  station  "P"  for  optical  parameterization  schemes.  Scheme  1  -  solid;  Scheme  2  - 
dashed;  Scheme  3  -  dotted. 
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Figure  13B.  Vertical  absorption  profiles  representative  of  waters  in  the  vicinity  of 
ocean  station  "P"  for  optical  parameterization  schemes.  Scheme  1  -  solid;  Scheme  2  - 
dashed;  Scheme  3  -  dotted. 


D.  APPLICATION  OF  OPTICAL  SCHEMES 

The  figures  presented  in  this  section  describe  the  results  of  applying  the  three 
optical  parameterization  schemes  during  different  seasons.  The  effects  of  the  individual 
schemes  are  compared  for  long  and  short  periods,  and  individual  deepening  and 
shallowing  events  are  examined  in  detail.  In  all  figures,  the  three  schemes  are 
represented  by  the  following  types  of  lines  ; 
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Scheme  1 
Scheme  2 
Scheme  3 


1 .  Short-Term  Seasonal  Comparisons 

Figures  14  through  17  illustrate  model  results  for  typical  10-day  periods  in 
each  of  the  four  seasons.  Part  A  of  each  figure  shows  the  atmospheric  forcing,  the 
resultant  uJj.,  and  mixed  layer  depth  during  the  period.  Part  B  depicts  the  evolution  of 
temperature  with  depth  throughout  the  period.  Part  C  shows  model  predictions  during 
a  subsection  of  the  lO-day  period  to  detail  specific  mixed  layer  changes  and  the  effects 
on  vertical  temperature  profiles.  The  color  contour  plots  of  part  D  portray  the 
temperature  difference  between  schemes  3  and  I  throughout  the  upper  50  meters  during 
the  period.  Scheme  3  is  assumed  to  be  the  most  realistic  of  the  three  schemes,  as  both 
the  spectral  attenuation  and  refraction  are  included.  Scheme  1  is  assumed  to  be  the  least 
accurate. 

The  differences  in  uJ?-  among  the  three  cases  is  most  apparent  when  Iq  is 
maximum.  The  greater  absorption  near  the  surface  for  schemes  2  and  3  causes  a  larger 
value  of  uJ>.,  which  leads  to  a  reduction  in  the  amount  of  turbulent  kinetic  energy 
available  for  entrainment.  The  effect  of  this  reduced  turbulent  entrainment  is  most 
evident  in  the  Figure  16A  portrayal  of  mixed  layer  depth,  as  a  smaller  value  of  h  is 
produced  by  shallowing  events  for  schemes  2  and  3.  For  mixed  layer  deepening  events. 
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the  decreased  entrainment  caused  by  a  greater  uJ}.  is  often  overcome  by  the  reduced 
thermal  gradient  being  entrained  for  schemes  2  and  3.  This  process  is  evident  in  Figure 
16A  where  scheme  2  displays  the  most  rapid  deepening  due  to  a  smaller  thermal  gradient. 
The  smaller  thermal  gradient  below  the  mixed  layer  results  from  the  increased  deep 
absorption  of  solar  radiation  and  decreased  intermediate  absorption  discussed  previously 
(Figure  13B). 

The  temperature  evolution  over  these  periods  also  reflects  the  effects  of  the 
absorption  curve  differences  for  the  three  schemes.  When  h  is  small,  the  increased  near 
surface  absorption  for  schemes  2  and  3  generates  greater  values  of  sea  surface 
temperature  (SST),  but  for  deeper  values  of  h  the  decreased  absorption  at  intermediate 
depths  dominates  and  SST  is  reduced  for  schemes  2  and  3.  This  effect  is  most  evident 
in  Figure  17C. 

The  increased  deep  absorption  for  schemes  2  and  3  produces  increased 
warming  at  depths  below  the  penetration  of  turbulent  entrainment.  At  50  meters  in 
Figure  15B  this  effect  over  the  10-day  period  produces  a  temperature  difference  of 
approximately  0.05°  C  between  .schemes  1  and  2.  At  these  depths,  solar  radiation 
produces  the  only  temperature  change  in  this  one-dimensional  representation.  This 
suggests  that  over  a  seasonal  or  annual  period  a  greater  amount  of  heat  is  imparted  to  the 
ocean  below  the  mixed  layer  than  represented  by  scheme  1. 

Significant  differences  in  temperature  at  specific  depths  develop  as  the  mixed 
layer  deepening  and  shallowing  events  occur  at  different  intervals  throughout  the  peritxl. 
The  jump  in  temperature  for  scheme  1  at  20  meters  in  Figure  15B  precedes  the 
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corresponding  increase  for  schemes  2  and  3  by  approximately  7  hours.  The  profiles  of 
Figure  16C  illustrate  that  the  differences  in  thermocline  gradient  for  the  three  schemes 
can  result  in  significant  temperature  differences  at  depths  within  the  thermoclines. 


The  color  contour  plots  of  temperature  difference  display  of  all  the  above 
phenomena.  The  differences  difference  between  schemes  3  and  1  are  evident  in  all  of 
these  figures.  Also,  the  dominant  effect  of  the  negative  absorption  difference  at 
intermediate  depths  over  the  smaller  positive  difference  at  shallower  depths  can  be 
observed  when  the  mixed  layer  significantly  penetrates  into  the  intermediate  depths.  The 
largest  contour  values  represent  the  temperature  gradient  differences  that  occur  in  the 
thermocline  regions  as  described  above. 

Overall,  it  can  be  seen  that  the  short-term  effect  of  the  different  schemes  is 
significant  during  all  seasons;  however,  the  predominant  methods  through  which  the 
differences  are  effected  change  due  to  the  change  in  seasonal  mixed  layer  evolution.  The 
first  few  days  after  model  initialization,  the  differences  between  schemes  are  relatively 
small,  suggesting  that  the  accuracy  of  optical  characterization  may  not  be  of  primary 
importance  to  short  range  predictions  of  the  acoustic  environment  based  on  thermal 
structure. 

2.  Seasonal  Cycles 

Figures  18  and  19  represent  the  evolution  of  the  thermal  structure  for  the 
different  schemes  during  seasonal  warming  and  cooling  cycles,  respectively.  The 
different  schemes  affect  the  seasonal  thermal  structure  by  the  same  mechanisms  as  for 


52 


shorter  periods.  The  long-term  results  are  due  to  a  succession  of  short-term  differences. 
The  profiles  in  Figures  18B  and  19B  demonstrate  the  typical  evolution  of  the  seasonal 
thermocline  during  these  periods,  as  well  as  the  difference  in  this  evolution  as  a  result  of 
the  optical  parameterization  scheme  used.  The  contour  values  of  Figures  18C  and  19C 
display  the  significant  differences  in  temperature  throughout  the  period. 

3.  Consideration  of  Realistic  Atmospheric  Conditions 

Figure  20  shows  the  results  of  including  cloud  cover  in  model  forcing.  The 
lesser  values  of  h  throughout  the  period  reflect  the  overcast  conditions  that  were  observed 
most  of  the  period.  This  analysis  is  used  to  demonstrate  that  significant  differences 
between  the  schemes  exist  when  the  long  wave  radiation  budget  of  actual  cloud  conditions 
is  included  in  model  forcing  parameters.  Schemes  2  and  3,  in  their  current  form,  cannot 
be  construed  as  accurate  for  these  conditions  as  realistic  effects  of  cloud  droplets, 
precipitation  and  water  vapor  are  not  parameterized. 

The  results  indicate  that  the  difference  between  the  schemes  for  this  10-day 
period  remains  perceptable  with  a  more  realistic  long  wave  heat  budget  included  in  model 
forcing.  Similar  to  the  previous  clear  sky  analyses,  more  absorption  occurs  in  the 
shallow  and  deep  regions  for  schemes  2  and  3,  and  less  absorption  occurs  at  intermediate 
depths.  This  causes  the  Scheme  2  and  3  predictions  to  reflect  warmer  SST  when  h  is 
small,  and  a  cooler  SST  when  h  is  large.  The  increased  warming  for  schemes  2  and  3 
at  60  meters  reflects  the  greater  absorption  at  depth. 
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Figure  14B.  Temperature  response  time  series  at  various  depths  during  a  typical  10- 
day  peri(xJ  in  Spring. 
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Figure  14C.  Close-up  of  mixcxl  layer  lemperature  and  depth  evolution  during  days 
6  and  7  of  the  Spring  period. 
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Figure  !4D.  Color  contours  of  temperature  difference  between  scheme  3  and  scheme 
I  for  a  typical  10-day  period  in  Spring 
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Time  (days) 


Figure  I5A.  Time  series  of  atmospheric  forcing  and  mixed  layer  response  for  a 
typical  10-day  period  in  Summer. 
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Figure  15B.  Temperature  response  time  series  at  various  depths  during  a  typical  10- 
day  peri(xl  in  Summer. 
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Figure  15C.  Close-up  of  mixed  layer  temperature  and  depth  evolution  during  day  8 
of  the  Summer  period. 
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Figure  15D.  Color  contours  of  temperature  difference  between  scheme  3  and  scheme 
1  for  a  typical  10-day  period  in  Summer 
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Figure  16A.  Time  series  of  atmospheric  forcing  and  mixed  layer  response  for  a 
typical  10-day  period  in  Fall. 
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Figure  16C.  Close-up  of  mixed  layer  temperature  and  depth  evolution  during  days 
5  and  6  of  the  Fall  period. 
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Figure  16D.  Color  contours  of  temperature  difference  between  scheme  3  and  scheme 
I  for  a  typical  10-day  pericxl  in  Fall  ). 
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Figure  17A.  Time  series  of  atmospheric  forcing  and  mixed  layer  response  for  a 
typical  10-day  period  in  Wintci. 
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Figure  17B.  Temperature  response  time  series  at  various  depths  during  a  typical  10- 
day  period  in  Winter. 
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Figure  17C.  Ciose-up  of  mixed  layer  temperature  and  depth  evolution  during  days 
7  and  8  of  the  Winter  period. 
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Figure  17D.  Color  contours  of  temperature  difference  between  scheme  3  and  scheme 
!  for  a  typical  lO-day  period  in  Winter  ). 
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Figure  18B.  Vertical  temperature  profiles  near  the  beginning  and  end  of  a  typical  60- 
du>  period  during  the  warming  season. 
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Figure  18C.  Color  contours  of  temperature  difference  between  scheme  3  and  scheme 
1  for  a  typical  60-day  period  during  the  warming  season 
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Figure  19A.  Temperature  response  time  series  at  various  depths  for  a  typical  60-day 
period  during  the  cooling  season. 
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Figure  19B.  Vertical  temperature  profiles  near  the  beginning  and  end  of  a  typical  60- 
day  period  during  the  cooling  season. 
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Figure  19C.  Color  contours  of  temperature  difference  between  scheme  3  and  scheme 
1  for  a  typical  60-day  period  during  the  ctxding  season  i  ) 
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typical  lO-day  pericxl  in  Fall  with  clouds  included  in  atmospheric  forcing. 
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Figure  20B.  Temperature  response  time  series  at  various  depths  during  a  typical  10- 
day  period  in  Fall  with  clouds  included  in  atmospheric  forcing.. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  purpose  of  this  study  was  to  more  accurately  prescribe  the  effect  of  solar 
radiation  on  mtxlels  for  the  thermal  structure  of  the  upper  ocean.  While  it  is,  as  yet, 
impossible  to  verify  the  degree  of  improvement  in  the  accuracy  of  predictions,  the  new 
parameterizations  do  make  a  difference  in  the  model  predictions  of  the  thermal  structure. 

The  nature  of  the  optical  interactions  between  sunlight  and  the  ocean  indicates  that 
both  spectral  attenuation  and  refraction  of  direct  sunlight  are  phenomena  that  should  be 
included  in  mixed  layer  modeling.  Optical  parameterization  scheme  3,  employing  the  use 
of  1 1  discrete  spectral  bands  and  considering  refraction  of  direct  sunlight,  should  be  the 
most  accurate,  provided  that  accurate  values  of  a,  and  6,  are  obtained  for  the  appropriate 
atmospheric  conditions.  (The  value  of  a,  prescribes  the  amount  of  total  insolation  in  each 
spectral  band,  and  R,  prescribes  the  amount  of  each  spectral  band  that  is  direct  sunlight.) 
This  scheme  requires  no  modification  to  use  spectral  attenuation  coefficients  for  specific 
liKations  measured  U)cally  or  remotely  by  satellite  colorimeters. 

The  numerical  case  studies  indicate  that  significant  changes  develop  in  the  thermal 
structure  because  of  the  complexity  of  optical  factors.  It  is  interesting  that  the  different 
optical  schemes  dt)  not  cause  consistent  changes  in  the  thermal  structure.  The  mixed 
layer  depth  is  not  always  more  shallow  or  more  deep;  SST  is  not  always  warmer  or 
c(H)ler.  The  complexity  of  the  overall  effect  is  attributable  to  the  inherent  complexity  of 
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the  relationship  between  thermal  and  turbulent  interactions  in  dynamics  of  the  ocean 
mixed  layer. 


The  differences  demonstrated  in  the  numerical  simulations  justify  the  consideration 
of  detailed  optical  otfecLs  in  regional  or  global  heat  budget  predictive  models.  The 
demonstrated  differences  in  temperature  distribution  over  seasonal  time  scales  could  be  * 

important  in  defining  the  (Kean's  role  as  a  heat  sink  or  source  when  coupled  to  a  dynamic 
atmosphere.  However,  evaluation  of  the  short-term  numerical  results  indicate  alternate 
optical  narameterizations  are  not  as  significant  for  predictions  of  the  thermal  structure 
to  be  applied  to  acoustic  forecasts.  The  thermal  .structure  projections  for  72-hour  periixis, 
typical  for  acoustic  forecasts,  do  not  have  differences  large  enough  to  significantly  alter 
acousiic  conditions,  especially  considering  the  scale  of  other  inaccuracies  involved  in 
using  a  one  dimensional  model  to  represent  acoustic  conditions. 

The  improved  optical  parameterizations  developed  in  this  study  repre.sent  a  small 
but  important  step  in  improving  the  accuracy  of  diagnostic  and  prognostic  mcxlels  for 
(Kean  mixed  layer  processes  through  more  accurate  repre.sentation  of  physical  prcKesses 
involved. 

B.  RECOMMENDATIONS 

The  following  steps  are  recommended  as  in  the  further  development  of  optical 
parameterization  schemes;  c 

•  Improve  the  parameterizations  of  atmospheric  effects  on  solar  radiation 
incident  at  the  (Kean  surface.  The  methods  described  in  chapter  II  to  (.alculate  /o  make 
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many  assumptions  regarding  the  effect  of  atmospheric  constituents.  Evaluation  of  the 
optical  parameterization  schemes  developed  in  this  study  use  similar  as.sumptions.  The 
use  of  models  describing  the  effect  of  these  atmospheric  constituents  on  solar  radiation 
should  be  included  to  more  accurately  specify  the  irradiance  at  individual  wavelengths. 
Satellite  observations  and  regional  or  global  data  fields  provide  sufficient  descriptions  of 
the  atmosphere  to  incorporate  these  models  into  irradiance  calculations.  Included  in  this 
is  the  need  to  incorporate  predictions  of  atmospheric  transmittance  used  in  the 
determination  of  albedo. 

•  Perform  verification  of  the  accuracy  of  these  improved  optical 
parameterization  schemes  at  specific  locations.  After  atmospheric  effects  are  more 
accurately  specified,  the  predictions  of  mixed  layer  models  that  incorptirate  improved 
optical  considerations  should  be  compared  with  the  same  mixlels  using  more  primitive 
parameter izations.  This  would  require  reevaluation  of  empirical  constants  in  the  models 
which  are  often  adjusted  to  provide  more  accurate  predictions  at  a  specific  location. 

•  Consider  the  vertical  and  tempiiral  variation  in  ocean  optical  properties. 
As  discus.sed  in  chapter  II,  proces.ses  such  as  phytoplankton  productivity  cause 
considerable  vertical  fluctuations  of  the  attenuation  characteristics  of  the  ocean  with 
season.  This  variability  may  result  in  layers  of  increa.sed  absorption  that  significantly 
effect  the  vertical  thermal  structure.  Resultant  short-term  variation  in  SST  may  be 
important  in  quantifying  the  air-sea  heat  exchanges  that  trigger  major  atmospheric  and 
oceanic  synoptic  scale  events. 
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